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SUMMARY 


Coastal upwelling is examined as it relates to the cycling of chemical species 
in coastal waters along the west coast of the North American continent. The 
temporal and spatial features of upwelling phenomena in the Eastern boundary regions 
of the North Pacific Ocean are presented and discussed in terms of upwelling epi- 
sodes. Climate conditions that affect upwelling include: thermal effects, wind- 

induced shear stress which moves surface layers of water, and the curl of the wind 
stress vector which is thought to affect the extent and nature of upwelling and the 
formation of offshore convergent downwelling fronts. These effects and the interac- 
tion of sunlight and upwelled nutrients which result in a biological bloom in sur- 
face waters is modeled analytically. The roles of biological and chemical species, 
including the effects of predation, are discussed in that context, and relevant 
remote sensing and in situ observations are presented. Climatological, oceano- 
graphic, biological, physical, chemical events, and processes that pertain to bio- 
geochemical cycling are presented and described by a set of partial differential 
equations. Simple preliminary results are obtained and are compared with data. 

Thus a fairly general framework has been laid where the many facets of biogeochemi- 
cal cycling in coastal upwelled waters can be examined in their relationship to one 
another, and to the whole, to whatever level of detail or approximation is warranted 
or desired. 


SYMBOLS 

a coefficient in eq. (12), defined by eq. (19) 

B biomass concentration 

B threshold biomass 

b coefficient in eq. (13) 

C D drag coefficient 

c mass fraction of species i, proportionality coefficient for wind velocity, 

eq. (16) 

Cp specific heat at constant pressure 

D diffusion coefficient 
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KV 


d particle size 

d e depth of euphotic layer 

f Coriolis parameter 

g gravitational constant 

Ip intensity of radiation backscattered upward (-z direction) 

If intensity of radiation transmitted in +z direction in sea 
i small distance from coast where properties are known 

J| negative rate of predation of species i per unit mass 
K v Kubelka-Munk absorption coefficient at frequency v 
k eddy conductivity 

unit vector directed vertically 

D 

P m specific biomass production rate at optimal light intensity 

fi Ekman mass transport, eq. (10) 

m slope of linear variation of wind, eq. (16) 

n index of refraction, compatibility condition fraction 

p pressure 

Q r defined by eq. (25) 

q r radiative flux 

R component of frictional force per unit volume 

r m maximum ration for zooplankton grazing 

Sj scattering coefficient at frequency v 

t time 

U air velocity in x-direction 

u mean horizontal velocity in water, eq. (35) 
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u water velocity component in x-direction 

V air velocity in y-direction 

$ horizontal air velocity (fig. 8) 

^adj adjusted horizontal air velocity (fig. 8) 
v water velocity in y-direction 

v g water velocity in eastward direction 

w water velocity in z-direction 

x,y orthogonal sea surface coordinate directions 

z . coordinate measured downward from sea surface 

a initial slope of specific photosynthesis rate at irradiance near zero 

e distance offshore where weather data are known 

6 distance offshore where convergent downwelling front is located 

$ latitude 

u radiative absorption coefficient 

p density, reflectivity 

0 zooplankton grazing rate constant 

t shear stress vector 

o production rate of species i, per unit volume 
Q angular velocity of rotation of Earth 

Subscripts: 
a air 

B blackbody 

H horizontal component 

1 species i 
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o value at, or near, coast, value at ocean surface 

pm particulate material 

r radiative, or residence time 

s surface value, scattering 

t terminal (settling velocity, eq. (20)) 

v vertical component 

x component in x-direction 

y component in y-direction 

z component in z-direction 

v spectral frequency 

6 value at offshore front 


INTRODUCTION 


Scientists tend to see things in bits and pieces, and understandably so. 
Creation is so complicated and is interdependent in such subtle Mays that it is a 
considerable task to become proficient in one narrow aspect of scientific study. 
This may be especially true in the marine sciences. With respect to biogeochemical 
cycling in Eastern boundary currents, or coastal upwelling, it is no longer suffi- 
cient to look at bits and pieces unrelated to the whole that is there. Tbe effects 
of climatology, oceanography, biology, physical and chemical events, and processes 
should be examined together. It would be useful to develop a framework which 
relates all these parts to see where they fit, in order to begin to assess their 
relative importance (some may not be important at all), and to ask or answer ques- 
tions concerning events that affect the entire ecosystem. 

Although runoff of water from landmasses into the oceans is a contributing 
factor toward the overall chemical budget in the seas, coastal upwelling is gener- 
ally considered to be the dominant mechanism for biogeochemical cycling along much 
of the west coast of North America. In the upwelling process, cold nutrient-rich 
water is drawn from the depths and is exposed to sunlight in the euphotic zone. 
Phytoplankton seed germinates there and grows by uptake of upwelled nutrients 
(derived in part from bacterial decomposition of organic matter on the shallow 
continental shelf), and it in turn is grazed upon by zooplankton in a complicated 
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interrelated process. The process is highly variable both temporally and 
spatially. 

Nevertheless, extensive literature exists about the various facets of the 
problem and we can formulate our thinking in specific and useful ways. In the 
context of this, study, many disciplines are involved: fluid dynamics; radiative 

transfer; transport of, and finite rate changes in, chemical specie concentrations; 
photo and biochemical processes; satellite and airborne sensing; and a host of 
others. Clearly, we must avoid trying to describe any one discipline in such detail 
that precludes a useful description of the overall biogeochemical system. 

For present purposes, we make every reasonable simplification to arrive at a 
description of the gross behavior of the system, and to reveal particular gaps in 
our understanding where useful research can be defined that would advance our under- 
standing of biogeochemical cycling in the sea. We begin by considering the spatial 
and temporal gross features of upwelling. 


COASTAL UPWELLING ALONG THE WEST COAST OF NORTH AMERICA 


Ocean and weather phenomena are interdependent. Figure 1 shows a modified 
National Oceanographic and Atmospheric Administration sea surface thermal analy- 
sis. The solid lines show temperature contours for 1984. Part a of the figure 
corresponds to early summer (May 8), at which time a weak thermal front (considered 
to be upwelled water) existed along the coast from Point Conception to Point 
Arena. Part b corresponds to winter (December 7). The sea surface temperature is 
lower, and there is no thermal front noted. We have added to that figure (obtained 
from ref. 1) to show isotherms from comparable times of the year 9 yr earlier 
(interpolated from ref. 2), as shown by the dashed lines. Generally, the ocean 
appears to be 1° to 2°C warmer in 1984 than it was in 1975. This may be an artifact 
from the warm El Nino current which first appeared in May of 1982, and was still 
evident in 1983 and part of 1984. That current raised the ocean temperature as much 
as 5°C, and globally wreaked havoc with the weather. The amount of energy repre- 
sented in 1° temperature change is enormous; the equivalent of raising the body of 
water to an altitude of 427 m. . 

As indicated by figure la, the dominant upwelling on the west coast of North 
America occurs in the spring and summer of the year. Long-term climatic data show 
that the wind tends to have a component that is directed toward the equator during 
that time. The wind, combined with the Coriolis effect of the rotation of the Earth 
tends to move surface waters in an offshore direction. These waters are replaced by 
cooler water drawn from the depths, which brings the ingredients for a biological 
bloom into the sunlit surface layers. At any latitude, the upwelling is intermit- 
tent as the wind varies, but the general effect is one which produces a biomass that 
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(a) May 8, 1984 solid contours; May 1 to 15, 1975 (dashed contours interpolated, 

ref. 2). 


Figure 1.- Sea surface thermal analysis. 


6 



SFOKAR «J 



DEC 1-5, 1975 DATA (INTERPOLATED NOAA/NMFS, REF. 2) 


(b) Dec. 7, 1984 solid contours; Dec. 1 to 15, 1975 (dashed contours interpolated, 

ref. 2). 

Figure 1.- Concluded. 
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moves sluggishly offshore. Figure 2 is a computer-enhanced infrared satellite 
photograph of upwelling along the coast of Northern California (ref. 3). The light 
tones indicate cool upwelled water, which forms irregular patterns or tongues and 
swirls as it moves offshore and mixes with the coastal California current. 



Figure 2.- Computer-enhanced infrared satellite image of California coast, May 5, 
1984 (NOAA). Light tones along the coast are cool, upwelled water. 


Biologically, the coastal upwelling is very important. A record marketable 
world fish catch of 80 million tons was made in 1984 (ref. 4). It has been esti- 
mated that half of the world's fish catch is obtained from 0.1 t of the ocean area — 
that which is affected by coastal upwelling (ref. 5). 
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The intensity of upwelling varies both spatially and temporally. Figure 3 
shows the 20-yr mean monthly upwelling index (defined subsequently) as a function of 
time for two coastal locations; the solid line corresponds to Point Sal (refs. 6, 7) 
which is near Point Conception; and the dashed line corresponds to near Cape Blanco 
in Southern Oregon (ref. 6). Both locations exhibit a peak in late spring or in 
early s umm er, and a minimum near the beginning of the calendar year. However, the 
long-term mean upwelling at Cape Blanco occurs only from about March to October — the 
ocean actually downwells at the coast during the balance of the year. By contrast, 
the upwelling at point Sal is more intense (the peak by a factor of 2.5) and is more 
persistent — the long-term mean shows dominant upwelling throughout the year. At 
this same location, the mean hourly upwelling index computed for April 4 to 13, 1981 
(ref. 7) is shown by the solid round symbol. 



Figure 3.- Upwelling index comparison. 

The detailed hourly variation of the upwelling wind (assumed to be the equator- 
directed component that is locally parallel to the coast (ref. 7)) is shown by the 
Solid line in figure 4. The windspeed is highly variable, with peaks near 12 m/sec, 
and with the minimum near zero on April 11 (the direction was actually reversed — 
which was briefly conducive to downwelling) . At the sea's surface, wind varies not 
hourly, but continuously. The variable upwelling wind causes the sea's surface 
temperature to rise and fall in response (dotted line, fig. 4), trending downward 
until the upwelling episode weakens (about April 13) as the wind level diminishes. 
The general decline in temperature during the episode is indicative of cooler water 
being drawn from the depths into the surface layers of the sea. 
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WIND SPEED COMPONENT FROM 320 deg, m/sec 

- SEA TEMPERATURE, °C 



Figure 4.- Point Sal buoy data (NOAA). 


The result of many such episodes is a large upwelling structure that persists 
far out to sea — well beyond the continental shelf, as shown by the light tones in 
figure 5. Some biological observations made in the "sample area" shown on the 
figure will be discussed subsequently. 

The variation of local peak of the long-term monthly upwelling index is shown 
as a function of latitude along the west coast of North America by the solid line in 
figure 6. The peak occurs in April at the southern part of Baja California, Mexico; 
and progresses northward with time, as the upwelling season develops. From around 
Monterey Bay northward, the peak occurs in July. For the entire west coast of North 
America, the most intense and persistent upwelling on a long-term mean monthly basis 
occurs between San Diego and Point Conception. On a weekly basis (ref. 8), the 
maximum local peak for the year may move from 33°N to 39°N and shift from April to 
August in time. In 1968, 1970, and 1971 the weekly mean peaked for the year at 33°N 
in June, May, and June, respectively. For 1967 the yearly maximum local weekly 
upwelling index peaked at both 36°N and 39°N in the week of May 21. In 1969, 1972, 
and 1973 the weekly mean peak moved to 39°N (Point Arena — fig. 2) and occurred in 
August, April, and May, respectively. From the Oregon border northward, the long- 
term mean-upwelling index is significantly diminished, probably because of the 
prevailing offshore location of the atmospheric high-pressure system. Thus the Gulf 
of Alaska is not dominated by coastal upwelling but is dominated instead by the 
effects of the cyclonic gyre of the Northeast Pacific Ocean. 

The dashed line in figure 6 is the latitudinal variation of the long-term mean 
of the curl of the wind-stress vector (ref. 9) associated with the peak upwelling 
index line. If the curl of the wind-stress vector is positive, the component of the 
upwelling wind increases offshore, and the upwelled water tends to continue to move 
farther offshore (fig. 7). On the other hand, if the curl is negative the upwelling 
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Figure 5.- Computer-enhanced infrared satellite image near Point Sal, April 13, 
1981 (NOAA). Light tones along the coast are cool, upwelled water. 
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LONG TERM PREVAILING CONDITION 


SEA SURFACE WIND 



• WIND INCREASES OFFSHORE 

(OR POSITIVE CURL OF THE WINDSTRESS 
VECTOR) 

• UPWELLED WATER MOVES FAR OFFSHORE 
(200 km) 


TEMPORARY IMPORTANT VARIATION 


SEA SURFACE WIND 



• WIND DIMINSHES OFFSHORE 

(OR NEGATIVE CURL OF THE WINDSTRESS 
VECTOR) 

• PROMOTES LOCAL OFFSHORE "FRONTS" OF 
BIOLOGICALLY RICH WATER (50 km) 


Figure 7.- Effect of the curl of the wind-stress vector on the upwelling structure, 
(a) Prevailing offshore divergence of upwelled water, (b) Convergent front. 

wind diminishes offshore, and a convergent offshore upwelling front tends to form. 
That is, the upwelling is arrested at the front, downwelling may occur offshore, and 
(in the simple schematic) a closed cell of recirculated water may develop. 

Except near Point Eugenia, Mexico, and locally north of Vancouver Island, 
Canada, the long-term mean of the curl of the wind stress vector is positive in 
figure 6, which is conducive to the divergence of upwelled water offshore. It is 
notable that although the long-term mean values of the curl are small in the region 
between Point Conception and San Francisco and around the Oregon border, this curl 
is still positive. Thus the long-term mean condition is conducive to offshore 
divergence of upwelled water. The two conditions are illustrated in figure 7; 
figure 7a shows the prevailing offshore divergence of upwelled water, and figure 7b 
shows the variation which produces a convergent front. 

However, we must recall that upwelling occurs in episodes, and during an epi- 
sode the curl of the wind-stress vector is highly variable and can become negative, 
so that a convergent offshore front. could be formed at some stage of the episode. 
During the upwelling season, the climate of the west coast of North America is 
dominated by a high-pressure region offshore, where the wind tends to circulate 
around that high in a clockwise direction. This high pressure region varies in 
strength and is not stationary. During a blow, as the high pressure region moves 
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toward shore, the wind diminishes at sea first and subsequently may cause a negative 
curl to occur so that at the end of the blow an offshore front is formed. For exam- 
ple, depending on the location of the high, it may move toward the coast so that the 
wind direction increases from the north and then from the east. Thus, in the San 
Francisco Bay area, the "typical” springtime marine layer and fog pattern is inter- 
rupted and does not move as far inland. The local curl of the upwelling wind-stress 
vector may become negative temporarily, which promotes the formation of an offshore 
convergent front. Or, if the wind becomes easterly or southeasterly, the upwelling 
episode can be arrested entirely. As an observation, Howe (ref. 10) has left port 
in strong winds, but found that offshore the wind abated (negative curl of the wind 
stress vector). Near the Farralon Islands, the wind vanished and a convergent off- 
shore front was visible, marked by a line of floating marine organisms as far as the 
eye could see. Indeed, commercial fishermen have learned to use weather information 
to help locate such fronts. Thus, if the weather report says that there is strong 
wind at the San Francisco buoy, and light wind at the Farralon Islands, some people 
will fish the Farralons because there may be a temporary, but biologically rich, con- 
vergent front there. Moreover (ref. 10), on an occasion when strong upwelling in 
Monterey Bay was terminated by a reversal of the upwelling wind, fish buyers reported 
that their facilities were strained to capacity one day, and idle the next. 


THE EQUATIONS OF CHANGE 


It is convenient and useful to relate phenomena that are important to biogeo- 
chemical cycling in the sea in mathematical terms through the equations of change. 
These include the usual equations of motion of oceanography as a subset, as well as 
equations that describe solar illumination, and a description of chemical and bio- 
logical events to some level of approximation. We use a left-handed coordinate 
system x,y,z for convenience in the northern hemisphere with z directed down- 
ward. The corresponding water velocities are u, v, and w. The conservation of 
mass can be written as 


3o 3(pu) 3(pv) 3(pw2 . 0 

at ax ay az 


(D 


Similarly, the conservation of momentum in the three coordinate directions is 


Du au au au au _ „ 1 ap 1 „ 

— = — + u — + v — + w — = 2vQ sin d> - — rr- + — R 

Dt at ax ay az 3 v p ax p x 


( 2 ) 


Dv.lv av av av 

Dt at ax ay az 


-2uQ sin $ - 


1|£ + 1r 

p ay p y 


(3) 


Dw aw aw „ aw t , aw _ . 

— = — + u — + v — + w — = -2v Q sin 
Dt at ax ay az e 


. 112 . 

p az 


+ - R 
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The substantial derivative which appears on the left of equations (2), (3), and (4) 
contains both acceleration and inertial terms which are expressed as the central 
equality. On the right-hand equality, the first term is the Coriolis effect of 
deflection caused by the rotation of the Earth, the central term is the pressure- 
gradient effect, and the last term contains the viscous effects. The gravity term, 
g, appears in equation (4) only, which is a statement of the vertical component of 
forces. The energy equation is 


„ II 

pc p at 


aT aT aT 

puc — + pvc — + pwc — 
d ax d av d az 



(5) 


where the right side is the vertical energy flux divergence, k is the eddy conduc 
tivity, and for the case of absorption, emission, and scattering of radiation on a 
spectral basis, the divergence of the radiative flux is 


a_ 

az 
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In this expression, Iq. and Io correspond to the transmitted and internally 

1 v v 

backscattered radiative intensity of frequency v (ref. 11), which are described by 
a system of two flux equations of radiative transfer (which include emission terras, 
refs. 11-14). 
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(8) 


Climatalogical conditions provide the boundary conditions at z = 0 for equa- 
tions (2) through (8). 


Finally the conservation of each species is stated as 


ac. 
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ac. 

+ u— U 
ax 
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(9) 


The right side of equation (9) contains the chemical, photochemical, biochemical, 
and predation processes that relate to the formation or depletion of each specie, i. 

Boundary conditions for these equations will be discussed subsequently, where 
appropriate. In this form, the solution of these equations is not a particularly 
appealing task. So we now consider important simplifications for specific ocean 
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applications, and how these may be extended to the present problem. Consider first 
some simple concepts of ocean currents as derived from approximations of the momen- 
tum equations (2), (3), and (4). 


OCEAN CURRENTS 


Currents Caused by Mass Distribution 

Major ocean circulations are governed by the pileup of water in particular 
regions of the Earth. Correspondingly, the equations of motion are simplified by 
neglecting friction, the vertical component for the Coriolis term, transient and 
inertial terms (actually, the entire, substantial derivative (ref. 15); that is, all 
of the terms to the left of the second equality in equations (2), (3), and (4)). 

This leads to the hydrostatic equations, where horizontal velocities are related to 
ocean elevation differences and the rotation of the Earth. (Incidentally, the large 
gyre that dominates the ocean circulation in the Gulf of Alaska, mentioned in con- 
nection with figure 6, could be treated in this manner.) 


Wind Driven Ocean Currents and Upwelling 

Upwelling has been considered for most of this century to be the result of 
wind-induced shear stress on the ocean's surface combined with the Coriolis effect 
of the rotation of the Earth. It is commonly assumed (refs. 16-20) for wind-driven 
currents, that both the hydrostatic pile up of water and the acceleration terms 
(actually the entire substantial derivative) can be neglected in the momentum equa- 
tions (2) and (3). These reduce to a balance of the viscous and Coriolis terms (the 
first and last terms on the right sides of the equations). The viscous stresses are 
the product of the eddy viscosity and the gradient of the horizontal velocity with 
respect to depth. 

The resulting equations were originally integrated by Ekman (ref. 1 6 ) , and 
subsequently by many others. A discussion of primary features of the upwelling 
solution appears in reference 18. These solutions of the momentum equations include 
various approximations and assumptions regarding surface wind stress, shallow- or 
deep-ocean-bottom boundary conditions, and constant or variable eddy viscosity. 

Ekman 's solution for the deep ocean yields a spiral horizontal velocity variation in 
which the amplitude is a maximum at the surface and decreases with depth. The 
result is that surface water is transported at 45° to the right of the wind direc- 
tion in the northern hemisphere. Importantly, the average motion of the Ekman layer 
is at a right angle to the wind. This is the conceptual basis for much of the study 
of upwelling phenomena. 

The component of Ekman transport directed offshore, computed for an adjusted 
geostrophic wind (rotated 15° to the left and reduced 30 % by friction), has been 
termed the upwelling index (refs. 6 and 8). It is an estimate of the transport of 
surface water offshore. The Ekman mass transport is expressed as 
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fi = i ; x 5 = i » a c D |t|? . E 


( 10 ) 


The component directed offshore has the units volume per-second-per-unit length of 
shoreline. The term f is the Coriolis parameter (a function of latitude), t is 
the sea-surface shear-stress vector resulting from the adjusted geostrophic wind, 
and i< is a unit vector directed vertically (refs. 6 and 8). The right-hand equal- 
ity expresses the shear stress in terms of air density, p a , a drag coefficient, C D , 
and the adjusted geostrophic wind velocity, 

In consideration of the highly idealized simplifications of the momentum equa- 
tions, the various assumptions employed, and the very irregular boundaries of ocean 
currents and landmass in nature, we must realize that some of the notions we have 
about upwelling may be mathematical artifacts, and not physical reality. Although 
the physics retained in the momentum equations balance the Coriolis term with fric- 
tion (and other terms by some investigators) (e.g., ref. 21), the Coriolis effect is 
thought to be important for large-scale ocean circulations, rather than for small 
ones. Since coastal upwelling is a somewhat spotty phenomenon, it may be surprising 
that momentum solutions describe upwelling as well as they do. The upwelling index 
has been used extensively to study the gross features of the transport of coastal 
waters seaward, and the variation of that transport with latitude and time. 

However, for the balance of this paper, it is convenient to abandon this 
approach. Ekman upwelling described above was originally derived for the deep 
ocean; for a shallow ocean, Ekman found the surface current to be more aligned with 
the wind. The shallow ocean (along the continental shelf) and particularly the 
motion of the thin euphotic layer at the surface (fig. 8) is of biological inter- 
est. The coordinate system is aligned with x directed offshore, normal to the 
coast, y is directed locally parallel to the coast toward the equator, and z is 
depth measured from the sea's surface. A quasi-steady state is assumed, and any 
along-shore flow v (x,y,z).*is allowed, but 3v/3y is neglected. The continuity 
equation (1) becomes for incompressible flow 


3u 

3x 



(ID 


Near the shore the solution for upwelling and the offshore movement of water is 

u = ax , w = -az (12) 

respectively. Similarly, if a convergent front exists offshore at x = 6, the 
solution of equation (11) in that vicinity will be shown to be of the form 

u = b(6 - x) , w = bz (13) 
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For generality, wind speed and direction just offshore can be adjusted arbitrarily 
to account for the frictional boundary layer, and Coriolis or other effects; hence, 
the surface velocity and its components in the surface layer of the ocean where 
photosynthesis occurs, is taken to be proportional to the local adjusted wind veloc 
ity and its components (the right side of fig. 8). 


WEATHER DATA 



Figure 8.- Upwelling/wind schematic, bathymetry -37° 50' n. 


The adjusted wind is taken to be a linear function of the distance from shore; 

*adj * V 1 ♦ “ ) ( ' i ' 

Consider that real-time weather data exist near shore, and at a distance of e 
offshore. At a distance x > nc (where n is a fraction) the components of the 
wind velocity are 


U(x) = U Q (1 + mx) , V(x) = V Q (1 + mx) (15) 

where it may be noted that the curl of the component of the wind vector parallel to 
the coast is 



(16) 


The surface velocity of the water in the x-direction is 
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x > ne 


(17) 


u(x) = cU = cU Q (1 + mx) , 

That relation can be derived readily for the laminar flow of wind blowing over water 
by matching shear stress at the water/air interface between boundary layers in the 
air and water. The proportionally constant, c, is essentially the cube root of the 
ratio of the density viscosity product in the air, to that in the water. It varies 
from 0.025 to 0.03 for surface temperatures between 10°C and 20°C. For turbulent 
wind/ocean conditions, a simple proportionality can be argued; the constant 0.035 
has been used in the U.S. Geological Survey (USGS) oilspill trajectory model 
(ref. 22). 

At x = ne, the velocity given by equations (12) and (17) are equal and 

u(ne) = cUq( 1 + mne) = ane 
Thus the coefficient, a, is determined to be 

a = — = cUmfl + — ) 
ne 0 V mne' 


( 18 ) 


(19) 


where n will be determined subsequently. Thus for a given wind condition, the 
motion of the surface waters and the curl of the wind vector can be approximated 
simply for present purposes. 

It may be noted in figure 8 that the horizontal scale is greatly fore- 
shortened. The bathymetric profile shown crosses the Gulf of the Farallones at 
latitude ~37°50'N. The ocean floor is a broad area that slopes gently to the edge 
of the continental shelf. A significant feature is that the depth is not excessive 
(£100 m). In principle, bottom waters and nutrient rich sediments can be upwelled 
to the sunlit surface over much of the region. 


SOLAR RADIATION IN THE SEA 


The absorption of sunlight in the surface layers of the ocean is of primary 
importance to biogeochemical events. Most of the ocean is essentially a wet desert, 
with the exception of limited regions wherein nutrients are exposed to sunlight in 
surface layers where they participate in the production of a biomass. Thus regions 
of coastal upwelling are of major importance to the biogeochemistry of the Earth. 

Three features of the sea surface must be mentioned as background to what 
follows. First, although the upwelling may arise from depths of the order of 
100 m, the euphotic zone where sunlight interacts with organic matter is much more 
shallow — of the order of 10 m. It is estimated that 90% of all sunlight is absorbed 
in the upper 10 m of the sea (refs. 23 and 24). Coastal water absorption in the 
blue, green, yellow, and orange parts of the spectrum is two or three times that in 
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oceanic water, and infrared absorption occurs at the surface. Secchi disc observa- 
tions made by the author in 150 m depths off Point Sal, California, showed that the 
white Secchi disc could not be seen at depths greater than approximately 4 m. 
Importantly, light at a wavelength of 450 nm that is absorbed by chlorophyll pig- 
ments in phytoplanckton has an extinction coefficient of about 1 m , and is dimin- 
ished by 90£ at a depth of less than 2.5 m. 

Second, it should be noted that for an ordinary '•random” sea of 2 m or so, the 
superficial euphotic layer may be considered to be well-mixed. 

Third, the sinking rate of phytoplanckton is exceedingly small in the euphotic 
zone. This can be estimated by forming the ratio of the depth, z, where light with 
a wavelength of 450 nm is attenuated by a factor q_ /q_ to the settling velocity 

v v 0 

of a biomass particle of size d. The resulting residence time is 



where in this expression only, K v refers to the physical absorption coefficient and 
not to the Kubelka-Munk absorption coefficient (a factor of 2 different). For 
biologically rich coastal waters where K y is of the order of .1 m, the residence 
ime required for a particle to sink to a depth where light is attenuated by 90^, is 
of the order of 100 yr to 10 yr for small particles of 5 um- to 20 um-diam. For 
larger particles or clusters of 100-ym size, the residence time is half a year. For 
these particles, the bio mas s is diminished by predation, starvation, disease, old 
age, or ocean currents; but not by sinking. For large particles or clusters of 
millimeter size, the residence time in euphotic coastal waters is of the order of a 
day. In clearer oceanic water where light absorption is weak, the residence time is 
increased by a factor 1/K v (as much as the order of 10). 

Thus for assessing the effect of sunlight on the sea, we need only to consider 
a superficial layer in which phytoplanckton has a long residence time, which has a 
thickness of less than 10 m, and which is well-mixed. This is a useful and impor- 
tant simplification. With that in mind, let us turn to the energy equation (5), the 
equation for the radiative flux divergence (6), and the differential equations of 
radiative transfer in a scattering, absorbing, and emitting medium (eqs. (7) 
and (8)). 

For purposes of energy considerations, absorption is large when compared with 
scattering (K v >> S y ) for phytoplanckton. Thus the scattering coefficient, S v , and 
the backscattered intensity, I v , can be neglected (of course neither can be 
neglected for remote sensing signal analysis — a related problem that is readily 
treated). Radiant emission at 450 nm can be neglected, and the radiant intensity of 
photons of frequency, v, in the z direction at a depth z can be obtained from 
equation (7) and is 
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( 21 ) 
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where p_ is the surface reflection coefficient for photons of frequency, v, 
a v 

incident on the ocean. Correspondingly, the divergence of the incident flux q r 
at depth z is from equations (21) and (6) z v 
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If we sum the contributions over all frequencies, and average the radiation absorbed 
by integrating across the euphotic layer of thickness, d g , we obtain 




4 r K (1 - p ) 
v 0 v s v 



-K z 
v 
e 


dz 


(23) 


For the quasi-steady state in a well-mixed layer where temperature varies slowly in 
the z and y directions, the energy equation (5) combined with equation (23) yields 


puc q (1 - P )(e ^ e - 1) (24) 

P 3X d e r v Q V 

whose solution describes the temperature variation in the euphotic layer as a func- 
tion of the distance from the shoreline. Both the temperature variation and light 
intensity are important because these affect the rates at which reactions proceed. 
The growth rate of phytoplankton diminishes with decreased temperature (ref. 25), 
and the rate of metabolism of ocean organisms is approximately doubled by a 10°C 
increase in temperature (ref. 26). 

For convenience, define 


CL = 


r " pc d 

p e v=( 


-y d 

q r (1 - Pg )(1 - e v e ) 
V 0 


(25) 


and if i is a small distance from shore where the temperature Tq is known, 

x = ix , dx = i dx (26) 

Equation (24) can be written 



+ constant 


(27) 
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With reference to figure 8, equation (27) can be integrated to obtain 


( l - T <A i ln x 

V = a » ’ 


x < ne 


( 28 ) 


Similarly, from weather data at x = e 


/ T e - T \ . _L_ / Wmex 

\ Q / cU Q m * n VI + mxj * 


x > ne 


(29) 


At x = ne, the temperature is continuous so that from equations (28) and (29), 

c n m = tn (1_£_SS_) * (-JEL-) S£ (30) 

\ Q / 0 VI + mne/ M + mne/ 1 


This is the transcendental compatibility equation whose solution is n, the match 
point in surface velocity which is shown in figure 8 and mentioned in connection 
with equation (15). 

For the case of a convergent front offshore at x = 6, 6 can be taken to be 
the known weather data location (by extrapolation if necessary). Thus for conve- 
nience, 6 = e and m = -1/6 so that in the vicinity of the front 


cU 


u = eU 0 (l - f) = -“(«- x) = b (6 = x) 


(31) 


which defines b and is of the form of equation (13) as mentioned previously. 
However, as the front is approached, x ♦ 6, and T fi -► ®; this is an unpleasant 
mathematical artifact. So let us say that at x = (6 - i), T * T^ in equation (10) 
which affects the first term of the right-hand side of the compatibility of equa- 
tion (30). Thus for offshore convergence, the compatibility relation can be put in 
the form 

cU n 

— r- = in 6(1 - n) + r in n« (32) 

o n — i 

where 



6 = T- » 1 

1 


( 33 ) 
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Correspondingly, equation (27) can be integrated to yield 




X . 

1 n 


(34) 


For offshore convergence, the compatibility relation for n is plotted as a family 
of curves in figure 9, and is illustrated as follows. For a day of average cloudi- 
ness along the west coast in July (ref. 27), the left side of equation (32) can be 
estimated in a straightforward way. The summation, Q in equation (25) may be 
0.027 W/cm 2 averaged over a 24-hr day. (It may be higher than that during a blow 
when no cloud cover is present, and lower as the blow diminishes and the cloud cover 
increases . ) 




Figure 9.- Match point location (eq. (32)). 


Assume that d e = 10 m, - T Q is 4 K, U Q is 10 m/sec (a blow), and 6 is 
approximately 50 km. The left side becomes 4.33. If Tq is known at i = 1 km, 

6 is 50, and n is found to be 0.22 from figure 9. the corresponding temperature 
in the euphotic layer is shown in figure 10 as a function of distance offshore. 

It may be noted in figure 10 that the temperature rises sharply near shore 
where the surface velocity (eq. 12) is small, so that the residence time for heating 
is enhanced. Similarly, near the front (6 = 50 for this example), the surface 
velocity diminishes to zero (eq. (18)), and both residence time and temperature are 
enhanced. The inset in figure 10 is a N0AA sea surface thermal analysis plot for 
July 1984. It shows surface isotherms that appear to converge near the frontal 
location off Point Arena. As mentioned above in figure 6, such fronts have been 
observed by the author south of that region under these approximate conditions. The 
raw temperature data shown are from 10°C at x = 1 to 13°C behind the front. The 
14°C symbol at x/6 = 1.16 actually corresponds to the first isotherm west (left) 
of the front. 
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Figure 10.- Surface temperature/ (biomass) variation; 6 = 50, 6 = 50 km, n = 0.22. 


This straightforward treatment of the energy and radiative transfer equa- 
tions (5)-(8), combined with the continuity equation (1), provides usable results 
for sea surface layer considerations. A more elegant solution of the energy equa- 
tion (5), in which only the partial derivatives with respect to t and y were 
neglected, was obtained in the form of confluent hypergeometric functions by Howe 
(ref. 10). However, the result was more bother than it was worth. A further sim- 
plification by the author in which the second partial derivative with respect to 
z is neglected (the first derivative is retained) has also led to a solution of the 
energy equation. Although it is more tractable than that of reference 10, it is not 
yet particularly illuminating, and awaits further consideration. 

Before proceeding, it should be noted once again that our considerations are 
qualitative estimates; that climate and ocean phenomena vary from moment to moment, 
with latitude and shoreline topography; and that established upwelling patterns are 
altered by weather change and by the sluggish California current that moves irregu- 
larly offshore toward the equator. 
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BIOLOGICAL AND CHEMICAL CONSIDERATIONS 


In the overall biogeochemical cycle in the sea, chemical nutrients and sunlight 
produce plant blooms in surface waters which in turn support animal life from proto- 
zoa to the higher forms. Both plants and animals ultimately die, and they and 
animal excreta sink and are decomposed by bacteria to form minerals which are 
upwelled in sunlit waters during that season of the year, and thus continue the 
intermittent cycle. Some features of the overall process are shown schematically 
(highly simplified) in figure 11, which is an attempt to illustrate concepts 
described in references 26 and 28. Oceanic and boundary current waters with their 
biochemical content contribute to this simple picture. In this section, we will 
attempt to illustrate some aspects of this cycle — both by exhibiting some biomass 
observations and by obtaining an idealized, well-behaved solution to a biomass 
equation. 

Figure 5 showed the satellite image of upwelled water along the central 
California coast on April 13, 1981 — the result of strong northwesternly winds shown 
in figure 4. The latter figure showed the thermal response of the ocean to the wind 
episodes that gave rise to upwelling characterized by an upwelling index near the 
20-yr-mean-monthly peak near Point Sal (the solid symbol in fig. 3). 

In situ biological estimates were made by the author near Point Sal on 
April 13, 1981. A secchi disk and a correlation of NASA and NOAA data (Nolton, 
ref. 29) were used to estimate chlorophyll concentration. The results appear as the 
three solid square symbols in figure 12. Of the observations taken four times 
during the day, two were identical, so the three symbols represent four observa- 
tions. The concentrations obtained are at the high end of the data base for the 
correlation, approximately 10 mg/m^. 

Two to four weeks earlier, airborne remote sensing data (ref. 30) in the Point 
Sal region showed low chlorophyll concentrations — 2 % to 5 % of those shown above. 

Wind conditions obtained by the aircraft at an altitude of 500 ft varied greatly in^ 
strength and direction during that time, with southerly winds March 18 and 20 being 
conducive to onshore downwelling and followed by northwester lies March 22-24 and 
27-30. No biological interpretation was given. Clearly, extensive water sampling 
should be obtained at meaningful times and at locations in the upwelling system to 
enhance understanding of biological phenomena, and to calibrate both surface and 
remote sensing observations. 

With respect to predation on the surface biomass, daytime shipboard echo- 
sounder traces on station near Point Sal on April 13, 1981, show zooplankton near 
the ocean floor (trace B, fig. 13). The bottom of the ocean is shown by the dark 
outline at the lower part of the traces, above which the water column is shown from 
approximately 80 m to 150 m below the surface. Dark masses in the water column are 
schools of fish. Identification was by the captain of the vessel — a veteran fisher- 
man who was able to predict the content of the net (quantity and species) before the 
net was brought aboard. Near the bottom, gray areas are noted as shortbelly, 
shrimp, plankton, and krill; rock cod and hake are noted on the trace, and from 
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nh£, excreta, protein metabolism 

NHJ, AMINO NITROGEN 


Figure 11.- Biogeochemical cycling diagram of the upwelling water column. 
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Figure 12.- Chlorophyll-Secchi data and correlation, present Secchi measurements. 


approximately 100 m to 150 m, hake and sablefish are shown. In other portions of 
the record that are not shown, schools of young hake were found close to the bottom. 

The plankton shown are not the surface phytoplankton; rather, these are 
zooplankton, which consist of small crustaceans, Jellyfish, worms, mollusks, eggs, 
and other larval stages. These animals consist of floaters and weak swimmers, which 
are characterized by diurnal vertical motion. At night, such a trace would show 
them to be near the surface, supposedly grazing on the phytoplankton. In the 
morning, these small creatures return to the bottom and in turn are grazed upon by 
bottom dwellers. Around sunset, the zooplankton mass begins to rise at a rate 
approaching 7 m/min (ref. 3D, so that the migration is complete by 10 to 11 p.m. 
This appears to be the mechanism whereby bottom fish are nourished by a surface 
phytoplankton bloom caused by the upwelling of nutrients into surface layers where 
photosynthesis supplies the needed energy. 

The water column in the upwelling regions of the continental shelf (fig. 11) 
may be considered as a surface layer and a bottom layer of intense biological activ- 
ity. That is, the euphotic layer at the surface is dominated by the photosynthetic 
production of phytoplankton in the spring and summer. The bottom is dominated by 
intense bacterial activity in a thin layer of detrius where organic matter is con- 
verted into minerals. In between those two layers, a mixture of chemical and bio- 
logical phenomena proceed with far less intensity. Thus for convenience, we may 
perhaps spatially -separate the water column to capture the main feature of 
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Figure 13.- Typical echo-sounder traces in lower part of water column with species 

noted. 


biogeochemical cycling. Moreover, if the formation of chemical species near the 
bottom occurs to a large extent in the fall and winter, it may be possible to sepa- 
rate future considerations temporally as well. Thus concentrations of nutrients may 
accumulate near the bottom during part of the year, and may be upwelled subsequently 
to the surface, to participate in the phytoplankton bloom. The spatial and temporal 
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simplification of biological and chemical events Mould facilitate analysis of the 
system, but it remains to be seen whether or not such an approach is realistic. 

The microbiology of the bottom is exceedingly complex by itself. Indeed, it 
may not be described in rigorous detail in the near future. Perhaps it does not 
have to be (after millennia in the fire age, we still cannot show how carbon monox- 
ide is oxidized — but fire is necessary and useful to civilization). Bottom water 
sampling of an ocean region of interest over a period of time can probably tell us 
all that is necessary to understand and model the system. The same is true of mid- 
water and surface water sampling. 

With reference to figure 11, bacterial activity is apparently high near the 
ocean's surface, diminishes with depth, and then it increases enormously in the 
bottom ooze. A result of bacterial decomposition of matter is that ammonia, carbon 
dioxide, and hydrogen sulfide are formed. Nitrification or microbial oxidation of 
ammonia to form nitrite and nitrate is an exceedingly important process. The nitro- 
gen cycle in the sea is complicated and has multiple pathways. Nitrite and ammonium 
may be only transient intermediate species in the production of the thermodynamic- 
ally stable nitrate. The nitrate is the second most abundant nitrogen species 
overall, and with chlorophyll, the nitrate plays primary roles in surface productiv- 
ity. The peak nitrate concentration occurs below the euphotic zone where oxygen 
concentration is a minimum. The most abundant nitrogen species is N 2 which is near 
saturation levels but which is generally inert. Nitrogen also appears dissolved in 
the amino form NHl, in nitrous oxide, and in particulate nitrogen-bearing sub- 
stances. During the winter months the surface concentration of nitrate reaches a 
peak, which is exhausted during the upwelling season when surface blooms of phyto- 
plankton occur. Chlorophyll is the nitrogenous component of particulate material 
which absorbs light at 450 nm. Although chlorophyll concentration is generally a 
fraction of a microgram per liter in the open ocean, it may be as high as 40 yg/J, in 
biologically rich coastal upwelled water where it is of primary importance to pro- 
ductivity. Importantly, chlorophyll can be sensed remotely or in situ (by various 
methods), and nitrate can be obtained by in situ samples which can be analyzed in 
the laboratory with little deterioration because of its stability. 

Carbon dioxide is dissolved in the water from the atmosphere, and is produced 
by the respiration of all organisms. The bulk of the carbon may come from skeletal 
remains decomposed by bacterial action. Sulfur is biologically important; its 
concentration appears to be a maximum in the bottom ooze as part of the microbial 
decomposition process. It may be detected at low tide when the aroma of H 2 S is 
noticeable. 

The above are qualitative notions about the biochemical phenomena involved in 
coastal upwelling. All of the details of these events are not yet known, and it is 
not possible to generalize the treatment of the problem from one region to 
another. Nevertheless the purposes of this introductory paper are served in the 
sense that a framework is provided in which to study and integrate these many clima- 
tological, oceanographic, biological, physical, and chemical events, processes, and 
observations. 
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For present purposes, it may be of some interest to form an idea of the gross 
behavior of the upwelling system under idealized conditions, quasi-steady, not 
nutrient-limited, and with "smoothly varying" predation. Several attempts have been 
made to model the phytoplankton biomass mathematically as it grows in nourished 
sunlit surface layers and as it is diminished by predation. A model developed by 
Denman and Platt (ref. 32) is 
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which can be shown to be of the form of the species of equation (9). Here u and 
w are mean horizontal and vertical velocities, D H and D y are horizontal and verti- 
cal diffusion coefficients, the next to the last term is the species production rate 
as affected by light intensity I z (z), and the last term accounts for predation. 
Almost any mathematician who goes to sea, as the upwelling wind episode subsides 
enough to allow that venture, quickly loses interest in vertical derivatives, and 
second derivatives of any kind for first estimates in surface layers. Indeed, 
during the upwelling season, the swell is generally from a southerly direction (a 
result of storms in Antarctica), and the wind-driven sea is generally from a 
northerly direction. The consequences are a generally confused sea surface and 
considerable mixing in the upper layers. Thus for the quasi-steady state, and for 
locally well-mixed surface water where along-shore derivatives are negligible, 
equation (35) reduces to 


u SB = P m B tanh 
ax tn 



r( 1 - e 
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(36) 


where for convenience, B has been set to zero, and I z (z) has been replaced with 
I . The latter is the radiative intensity of frequency, v, which is considered to 
be completely absorbed in the well-mixed euphotic layer. It is thus the unreflected 
portion of the incident photons of frequency which takes part in photosynthesis. 

Further assumptions are made to derive a simple solution to equation (36) so 
that a baseline idea of how the upwelling biosystem behaves under favorable circum- 
stances can be obtained. These assumptions may, of course, be relaxed in a more 
detailed study. Thus for present purposes, if 
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Under these conditions, equation (36) becomes 


u ^ = BaT -roB 
dx v m 


(39) 


The solution yields the phytoplankton biomass variation as a function of distance 
offshore ; 


B(x) = B n e 
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dx/u 


(40) 


Again if the exponent of the exponential is sufficiently small, this becomes 


B(x) - B q 
B 0 (a! v " r m o) 



(41) 


which is the same form as equation' (27). (It may be noted that the integral in the 
exponential in eq. (40) is the time required to move surface waters offshore from 
x = i. Thus, for large x, the exponential may be sufficiently large that eq. (40) 
is more appropriate than eq. (41).) 

The solution of equation (41) for offshore divergence is identical to equa- 
tions (28) and (29) with the left side of equation (41) replacing those of equa- 
tions (28) and (29), appropriately. 

The solution of equation (41) for offshore convergence is the same as it is for 
equations (32) and (34) with similar replacements, and is shown in figure 10 with 
the left ordinate labeled as shown. 


The result for both the divergent and convergent cases is that the biomass 
concentration variation is identical to the surface temperature variation. Thus for 
this idealized baseline case, the biomass concentration variation could be sensed 
with a thermometer or estimated from NOAA sea surface thermal analysis biweekly 
charts . 

Of course, upwelling is not an idealized phenomenon. There are several factors 
which are implicit in equation (41): (1) sufficient nutrient availability is indi- 

cated; (2) spatial variation of predation is negligible; (3) steady state is 
assumed; (4) there is no intrusion of biologically poor oceanic water; and (5) there 
is an absence of other unpleasant influences. Biomass concentration variations are 
often irregular and are contrary to expectation for reasons which need to be 
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understood. But we have a notion of what upwelling could produce biologically, and 
a reinforced conviction that remote and in situ observations, laboratory studies, 
and further analysis which relate to these features are essential to monitoring and 
understanding the entire interrelated upwelling/biological system. 


CONCLUDING REMARKS 


Major spatial and temporal features of coastal upwelling along the west coast 
of North America have been presented and discussed. Some local biological and 
physical information obtained by buoy data and in situ observations which include 
the entire water column have been presented as well as large scale data obtained by 
satellite and aircraft remote sensing. Features of the overall biogeochemical 
cycling system as affected by coastal upwelling have been outlined and where possi- 
ble, have been related to the observations. Discrepancies have been cited and the 
need for information that will enhance the understanding of the processes that 
affect the cycling of chemical species in the waters along the continental shelf has 
been noted. Knowledge of both the features of the cycling phenomena and the rates 
at which they change needs to be enhanced. 

The relationships of weather, oceanography, biology, chemistry, light, and 
physics have been formulated deterministically. This formulation provides a conve- 
nient framework for performing research, a means for ranking the importance of 
specific biogeochemical features (by ordering and by comparing terms in the equa- 
tions of change), and as a means of identifying basic information that is needed to 
solve the cycling problem to any level of detail that is desired. As a first 
approximation, the main features of biogeochemical cycling in upwelling systems may 
be divided spatially into two parts for convenience: ocean surface phenomena which 

include the phytoplankton bloom, and ocean floor phenomena which include the micro- 
bial conversion of organic material to minerals. It may be possible to further 
divide the problem temporarily, the link in time being the onset of upwelling which 
brings bottom nutrients into the thin surface euphotic zone. 

A preliminary illustrative result for ocean surface layers has been obtained 
which includes the effects of weather, oceanography, dynamics, energetics, light, 
biomass production, and predation for a "well-behaved" biogeochemical cycling system 
in upwelling. The result has been compared with sea-surface tghermal -analysis 
observations near the California coast. A biomass-temperature analogy for this 
well-behaved system was shown. The analysis could be perturbed for irregular or 
"poorly behaved" effects. It was noted that the equations which describe light 
scattering and absorption in the sea can be solved quite readily relevant to remote 
sensing signals. Future research activities were identified. Weather data for 
selective locations should include information concerning the spatial and temporal 
variation of surface wind shear stress, the curl of the wind stress vector, and both 
incident and backscattered light intensity on a spectral basis. Oceanographically, 
selective surface water sampling and analysis should be obtained and related to the 
progression of upwelling events. Transient phenomena and offshore frontal formation 
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data should be acquired and related to the behavior of the system. Biologically, 
the spatial and temporal variation of bottom activity should be examined by selec- 
tive sampling of bottom water and sediments to sort out the dominant features of 
chemical and microbiological reactions and their rates as they are affected by the 
behavior of the overall upwelling system. With respect to remote sensing, spectral 
scattering in water by plankton and other particulates should be characterized. The 
effect of the atmosphere and the air/sea interaction on light scattering and attenu- 
ation should be examined more thoroughly. Corresponding solutions to the Kubelka- 
Munk equations of light transfer through densely distributed scattering sites should 
be obtained. Computationally, the modeling of biogeochemical phenomena should be 
improved. Transient and three-dimensional effects should be examined. A bottom 
flow analysis should be developed which includes the conversion of organic to inor- 
ganic species. 

In the long term, solutions to the entire set of the equations of change over 
the entire domain, including rich physical and biological detail, the proper rate 
and property inputs, and the appropriate boundary condition (including climatologi- 
cal variations) should be able to describe the state of the entire biogeochemical 
cycling system, and could be used predictively to estimate the future course of 
events. As techniques improve they have potential application to other upwell- 
ings: along the west coasts of Africa and South America, Polar ice shelves, and 

boundary currents of major offshore ocean circulations. It may be that the redeem- 
ing feature of biogeochemical cycling in the vast ocean is that dominant events tend 
to be rather localized. 
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